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Abstract 

An indentation testing method, which utilizes lateral contact of a long cylindrical 
indenter, is developed for a thin transversely isotropic incompressible elastic film 
deposited onto a smooth rigid substrate. It is assumed that the material symmetry 
plane is orthogonal to the substrate surface, and the film thickness is small compared 
to the cylinder indenter length. The presented testing methodology is based on 
a least squares best fit of the first-order asymptotic model to the depth-sensing 
indentation data for recovering three independent elastic moduli which characterize 
an incompressible transversely isotropic material. The case of a weakly compressible 
material, which is important for biological tissues, is also discussed. 
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1 Introduction 


In recent years, the atomic force microscopy (AFM) indentation has become a useful 
technique for quantifying the mechanical properties of soft polymers p5]J and biological 
materials [IS]. It is well known dD! that many biological tissues are characterized by 
anisotropy, and the case of transverse isotropy plays an important role for describing their 
symmetry properties m 

Recently, the so-called cylindrical lateral indentation test, which utilizes lateral contact 
of a cylindrical indenter (see Fig. [Tj) , was developed by [TBJ under the assumption that 
both Poisson’s ratios of the tested transversely isotropic elastic material are known. Such 
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a test can be approximately modeled [4] in the framework of the two-dimensional contact 
model for an orthotropic elastic strip, which was was studied previously in DEI- We note 
that in the asymptotic model [4] it is assumed that the him thickness h and the indenter 
radius, R , should be small compared to the cylinder indenter length, l, as well as the 
contact zone size, a, should be small compared to the elastic him thickness h. 

However, although it is usual to assume in the indentation testing of an isotropic material 
that its Poisson’s ratio is known in advance j6j, this assumption in the case of material 
anisotropy requires a careful consideration, especially for almost incompressible materials. 
On the other hand, in the case of incompressible transversely isotropic material, the 
previously developed model |4j needs to be further refined. 

The rest of the paper is organized as follows. In Section [2j the constitutive equations 
for a transversely isotropic material are written out in the principal coordinate system. 
The case of incompressible material is considered in Section [3j In Section [4| we outline 
the hrst-order asymptotic model for the displacement-force relationship. In Section |5j the 
material identification procedure for the cylindrical lateral indentation test is outlined 
and the corresponding anisotropic materials properties evaluation procedure is illustrated 
by an example. Finally, in Section [6j we present a discussion of the results obtained and 
formulate our conclusions. 


2 A transversely isotropic elastic film deposited onto a rigid substrate 


We assume that the tested tissue sample forms a homogeneous, transversely isotropic, 
and linearly elastic layer is deposited onto a smooth rigid substrate in such a way that 
the plane of isotropy is perpendicular to the substrate surface. If the aqa^-plane coincides 
with the plane of isotropy, Hooke’s law for a transversely isotropic material can be written 
as follows [T3j : 


£11 — ^(^ll _ z/cr 22) — ^ 0 ' 33 ’ £23 ~ 9C^°' 23 ’ 

1 , s V\\ 1 

£22 — U<J U + < 722 ) — £ 13 — 

V\\ , -I 1 

£33 — — — (cm + C22) + ^r 0 ^, £12 — 


( 1 ) 


Here, £ lJ are infinitesimal strains, cq-j are stresses, E and v are Young’s modulus and Pois¬ 
son’s ratio in the symmetry plane, G is the in-plane shear modulus, E» is the longitudinal 
Young’s modulus, v\\ is Poisson’s ratio for loading in the longitudinal direction, and Gm 
is the longitudinal shear modulus. Recall also that transversely isotropic materials are 
characterized by five independent elastic constants (namely, E, E\\, Gy, u, and v\\), so 
that we have G = E/[ 2(1 + v)\. 


We assume that indenter has a shape of long cylinder oriented parallel to the material 
surface and is made from rigid material, so that the indenter deformation is negligible. 
Because the layer thickness is assumed to be small compared to the cylinder indenter 





Figure 1. A homogeneous, transversely isotropic, and linearly elastic layer deposited onto a rigid 
substrate. Fibril direction is perpendicular to the cylinder indenter axis. 

length, we apply a two-dimensional approximation for the stress-strain state beneath the 
indenter. We consider two cases, which differ by the orientation of the indenter cylinder 
with respect to the longitudinal axis of the transversely isotropic material. Namely, in the 
longitudinal case the indenter cylinder is parallel to the axis of material symmetry, while 
in the transverse case it is parallel to the plane of isotropy and thus perpendicular to the 
material symmetry axis. 

If the indenter is oriented along the aq-axis (see Fig. [TJ, the strain-stress relations for the 
case of plane strain deformations in the ^ 2 ^ 3 -plane are obtained from Eqs. (JTJ) under the 
assumption that £u = £u = £13 = 0 in the following from: 
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If the indenter is oriented along the X 3 -axis, the strain-stress relations for the case of plane 
strain deformations in the ayny-plane are obtained from Eqs. ([Tj) under the assumption 
that £13 = £23 = £33 = 0 in the following from: 
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Equations (J3| and (j2]) correspond to the longitudinal (along to the a^-axis) and transverse 
(perpendicular to the ^ 3 -axis) orientations of the cylinder indenter, respectively. 


3 Tissue incompressibility assumption and fitting parameters 


If the tissue material is assumed to be incompressible, then the following condition must 
be satisfied: £n + £22 + £33 = 0. As was shown by [T2], this incompressibility condition 













imposes two additional constrains on the components of the compliance tensor. Thus, 
for an incompressible transversely isotropic material only 3 material constants remain 
independent and the following relations hold [HIE]: 


v \\ 


1 

2 ’ 



( 4 ) 


At the same time, the condition of positive definiteness of the compliance tensor reduces 
to 


E < AE\ 


( 5 ) 


Since Poisson’s ratio u\\ is known in advance, while Poisson’s ratio v is expressed in terms 
of Young’s moduli E and 2E \\, the problem of material parameters identification reduces 
to the evaluation of the following clastic constants: E, E», and Gy. However, instead of 
these three dimensional parameters, we will make use of the transverse Young’s modulus 
E and two dimensionless ratios 
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Thus after evaluating the primary fitting elastic parameters E, n, and k, the elastic moduli 
En and Gy can be easily determined by formulas 

E\\ = nE, G\\ = kE. (7) 


Note that in light of (|5|), we will have n > 1/4. 


4 Displacement-force relation in the cylindrical lateral indentation 


During the depth-sensing frictionless indentation, the indenter displacement w (indenta¬ 
tion depth) is related to the indentation load P (contact force) by some nonlinear relation, 
which is approximated by the following first-order asymptotic formula [4]: 
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Here, Wo is a characteristic length parameter, Pq is a characteristic force parameter. 


The characteristic length parameter Wo is determined solely by the geometry of both the 
elastic sample and the rigid indenter as follows: 
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Here, h is the film sample thickness, R is the radius of the indenter. 






The characteristic force parameter Pq is defined by the formula 
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where l is the length of the indenter, and depends on an elastic constant 6, which will be 
defined later, and in the case of indentation testing, it is not known a priori. 


Further, the asymptotic constants do and d\ are dimensionless parameters, which charac¬ 
terize not only the elastic properties of the film sample, but also the contact conditions 
between the sample and the rigid substrate. Under the assumption of frictionless contact 
between the film and the rigid substrate (i. e., u 2 = 0 and <t 12 = cr 32 = 0 at x 2 = —h), we 
can make use of the following analytical formulas BEI : 
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Here, L(u) is a dimensionless function, which depends on the elastic properties of the film 
sample as well as on the orientation of the indenter [lj. 


4-1 Case of transverse indenter orientation 


In the transverse case (indenter is oriented parallel to the aq-axis and perpendicular to 
the aq-axis), the asymptotic constants dg and df are evaluated according to equations 
(18) and (19), where the expression for the kernel function L(u) depends on the values of 
the elastic constants as follows E0): 
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I 11 the case of an incompressible material, the dimensionless parameters A 1 and B L are 
given by 

2 n 2 — (4 n — 1 )k 


kU = 


k(An — 1) 


B x = l, 


(16) 






















while the corresponding elastic constant 6 L is given by 
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We note that the notation A L and B L is consistent with the notation E» and reflects the 
transverse orientation of the indenter cylinder to the axis of material symmetry. 


4-2 Case of longitudinal indenter orientation 


In the longitudinal case (indenter is oriented along the X 3 -axis), we have = 1 and = 
1. Therefore, the asymptotic constants d [J and d}\ are independent of elastic properties of 
the sample and are evaluated as follows: 

4J = 0.3517, d\ = -0.521. (18) 

In this case for incompressible material, the elastic constant 6^ is given by formula 

2 nE 


9 " = 


3 ’ 


(19) 


where the dimensionless parameter n is defined by the first formula 


5 Material identification procedure 


In the transverse case, by fitting the experimental depth-load data with the analytical 
approximation Q in view of ( fllj ) and (12)—(15), one can determine the elastic constants 
n, m, and 

_ < 20 > 

where Pf- is the corresponding best-fit value for Pq in 
Thus, in view of ( Jl7| ) and (20), we obtain 
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Formulas (20) and (21) solve the formulated above material properties identification prob¬ 
lem. 


In the longitudinal case, because the asymptotic constants djj an d d)\ are independent 
of the material properties, fitting the experimental data for the force-displacement re¬ 
lationship with the analytical approximation (J8]) allows to determine only one elastic 
constant 6^. 











5.1 Example. Application to indentation testing of a collagen layer 


To illustrate the outlined identification procedure, we consider the same cylindrical lateral 
depth-sensing indentation test for a collagen coating freely deposited on a glass substrate 
with the fibril direction perpendicular to the cylinder axis (see Fig. [Tj) as considered in 
our previous study [4]. The collagen membrane thickness is approximately 5 /im, while 
the indenter cylinder length is about 18 /mi, and the cylinder diameter is 6 /mi. 


Under the assumption of material incompressibility, the following best-fit values were 
obtained: 6 1 - = 10.96 MPa, n = 5.18, k = 1.35. Correspondingly, the normalized root- 
mean-square error of the best-fit approximation is 5 = 0.021. Now, using the first two 
formulas (21), we get approximation for the elastic moduli and the longitudinal shear 
modulus as follows: E = 4.03 MPa, E\\ = 20.86 MPa, and G\\ = 5.4 MPa. In light of Q 
we have u» = 0.5 and v = 0.90. Note also that the dimensional parameter m = G\\/G, 
which is related to k through formula m = 2(1 + v)k, takes the following value: m = 5.13. 


Recall that under the assumption of prescribed values for the Poisson’s ratios u\\ = 0.45 
and v = 0.45, the following data was obtained 0: 9 1 = 10.60 MPa, n = 8.25, m = 1.05, 
E = 16.9 MPa, E\\ = 139.6 MPa, and Gy =6.1 MPa, whereas 5 = 0.020. Therefore, the 
complete incompressibility is a rather restrictive (strong) hypothesis, since the longitudinal 
elastic modulus E\\ demonstrated a 7 times decrease. Thus, in order to be sure that the 
estimated Young moduli are correctly evaluated, such an assumption should be thoroughly 
justified and a level of material incompressibility should be clarified. 


6 Discussion and conclusions 


As it follows from the previous discussion, one can expect that for weakly compressible 
materials the material identification data is sensitive to the choice of the Poisson’s ratio 
values. Taking for instance the typical values v\\ = 0.49 and v = 0.95 (assumed in [9j for a 
soft highly anisotropic tissue with E\\/E = 20.3) and utilizing the identification procedure 
based on the first-order model [4], we, in particular, obtain n = 1.57 and 6 = 0.329 (see 
Table [I]). Note that the best-fit approximation error turns out to be an order of magnitude 
higher than in the previous cases (see Section 5.1). Moreover, the obtained value of the 
elastic moduli ratio n — Eu/E comes in contradiction with the following thermodynamic 
limitation imposed on the Poisson’s ratio v m : 


1 < v < 1 
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( 22 ) 


We also present in Table [l] the result of evaluation of the Young moduli for different values 
of the Poisson’s ratio v = 0.7,0.75,0.8,0.85,0.90. The material parameters evaluated for 


v = 0.90 do not also satisfy the restriction (22) and thus cannot be accepted. 


Therefore, in the case of slight incompressibility similar to that considered in [9], it is sug¬ 
gested first to solve the material parameters identification problem using the algorithm 





Table 1 

Variation of the evaluated moduli ratio n = E^/E, the upper bound u* = 1 — 2z/j?/n for the 
Poisson’s ratio v, and the approximation error 5 for different values of the Poisson’s ratio v and 
a prescribed value of the Poisson’s ratio v» = 0.49. 


V 

0.95 

0.90 

0.85 

0.80 

0.75 

0.70 

n 

2.47 

1.57 

13.52 

7.54 

8.12 

8.21 

V* 

0.806 

0.694 

0.964 

0.936 

0.941 

0.942 

5 

0.742 

0.329 

0.139 

0.021 

0.020 

0.023 


proposed in this paper under the assumption of the complete incompressibility, thus ob¬ 
taining an upper bound for the Poisson’s ratio v. After that, utilizing the algorithm [4], 
one can adjust step-by-step the values of Poisson’s ratios in such a way that the right-hand 
side inequality in (22) is preserved, while the approximation error is kept to a minimum. 


It is also worth emphasizing that in the case of thin layers, when the size of the contact 
zone is much larger than the layer thickness, the concept of a weakly incompressible layer 
can be introduced |16j . which takes into account the interplay between the material and 
geometrical parameters. The same approach can be utilized in the indentation problem 
under consideration. However, in the case of anisotropic material, one should take care 
of two small parameters, which define how the two bulk moduli describe the material 
behavior when its incompressibility increases. 


Summarizing, in the present paper, the first-order asymptotic model for the frictionlcss 
unilateral contact of a cylindrical indenter is applied for the development of cylindrical 
lateral indentation tests for thin (compared to the cylinder indenter length) transversely 
isotropic elastic samples under the assumption of material incompressibility. The proposed 
indentation material parameters identification method utilizes the depth-load indentation 
data for the cylinder indenter oriented orthogonal to the fiber direction only. If the depth¬ 
load data are available for the other indenter position (along the fiber direction), they can 
be additionally used for verifying the predicted value of the transverse modulus E. 
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